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ABSTRACT 

This  work  utilizes  Improved  Solar  Observing  Optical  Network  continuum  (630.2  nm)  and  Ha  (656.2  nm)  data  to: 
(1)  detect  and  measure  intrinsic  sunspot  rotations  occurring  in  the  photosphere  and  chromosphere,  (2)  identify  and 
measure  chromospheric  filament  mass  motions,  and  (3)  assess  any  large-scale  photospheric  and  chromospheric  mass 
couplings.  Significant  results  from  2003  October  27-29,  using  the  techniques  of  Brown  et  al.,  indicate  significant 
counter-rotation  between  the  two  large  sunspots  in  NOAA  AR  10486  on  October  29,  as  well  as  discrete  filament 
mass  motions  in  NOAA  AR  10484  on  October  27  that  appear  to  be  associated  with  at  least  one  C-class  solar  flare. 

Key  words:  Sun:  chromosphere  -  Sun:  filaments,  prominences  -  Sun:  flares  -  Sun:  photosphere  -  sunspots 
Online-only  material:  animation 


1.  INTRODUCTION 

Solar  atmospheric  activity  has  often  been  suggested  as  a 
precursor  to  solar  flares,  which  can  assist  in  devising  techniques 
to  accurately  predict  space  weather  events.  Intrinsic  sunspot 
rotations  have  been  detected  for  decades  (Evershed  1910; 
St.  John  1913;  Abetti  1932;  Maltby  1964;  Gopasyuk  1965; 
Bhatnagar  1967;  McIntosh  1981;  Pevtsov  &  Sattarov  1985;  Ding 
et  al.  1987),  but  their  association  with  solar  flare  generation  is 
uncertain.  Filament  activity  prior  to,  and  during,  flaring  also 
occurs  frequently  (Jiang  &  Wang  2001;  Jing  et  al.  2004;  Yan 
et  al.  2012),  but  establishing  diagnostic  pre-flare  indicators 
remains  elusive. 

A  history  of  detecting  and  measuring  mass  motions  in  the  so¬ 
lar  atmosphere  exists,  which  includes  intrinsic  sunspot  rotations 
and  filament  motions  across  the  solar  disk.  Intrinsic  sunspot 
rotation,  defined  as  the  proper  rotation  of  a  sunspot  around  an 
axis  perpendicular  to  the  solar  surface  and  centered  on  the  um¬ 
bra,  has  been  detected  since  the  early  20th  century,  utilizing 
both  ground-  and  space-based  photospheric  and  magnetogram 
data  (Evershed  1910;  St.  John  1913;  Abetti  1932;  Maltby  1964; 
Gopasyuk  1965;  Bhatnagar  1967;  McIntosh  1981;  Pevtsov  & 
Sattarov  1985;  Brown  et  al.  2003;  Min  &  Chae  2009;  Yan  & 
Qu  2007;  Zhang  et  al.  2007;  Yan  et  al.  2008,  2009;  Kazachenko 
et  al.  2009,  2010). 

Recent  years  have  seen  additional  work  on  this  topic  thanks 
to  spacecraft  such  as  the  Solar  and  Heliospheric  Observatory 
( SOHO\  Domingo  et  al.  1995),  Transition  Region  and  Coronal 
Explorer  ( TRACE ;  Schrijver  et  al.  1996),  and  Hinode  (Kosugi 
et  al.  2007),  along  with  ground-based  observatories  such  as  the 
Improved  Solar  Observing  Optical  Network  (ISOON;  Neidig 
et  al.  2003). 

Despite  abundant  detections  of  rotating  sunspots  in  the  past 
century,  the  causes  and  implications  of  this  phenomenon  are 
uncertain.  Yan  et  al.  (2008)  conducted  a  statistical  study  of  182 
measured  sunspot  rotations  using  SOHO ,  TRACE ,  and  Hinode 
data  to  constrain  sunspot  rotational  properties.  While  Yan  et  al. 
(2008)  report  that  there  is  no  rotational  (clockwise  (CW) 


versus  counterclockwise  (CCW))  preference  in  the  northern 
and  southern  solar  hemispheres,  the  northern  hemisphere  has 
a  greater  abundance  (12%)  of  rotating  sunspots,  while  positive 
(negative)  polarity  in  the  northern  (southern)  hemisphere  is 
about  twice  as  common  for  rotating  sunspots. 

In  a  statistical  study  of  chromospheric  penumbral 
whorls,  as  traced  by  superpenumbral  filament  directions, 
Balasubramaniam  et  al.  (2004)  found  only  a  weakly  correlated 
sense  of  twist  of  penumbral  fields  for  following  hemispheric 
preference.  However,  their  data  were  not  intended  to  study  the 
temporal  evolution  of  these  twists  as  related  to  flares. 

Intrinsically  rotating  sunspots  appear  to  be  relatively  rare 
occurrences  as  only  ~5%  of  the  analyzed  active  regions  in  Yan 
et  al.  (2008)  exhibit  rotations.  The  occurrence  of  more  than 
one  rotating  sunspot  per  active  region  is  more  rare  (~1%),  with 
only  a  minor  preference  for  multiple  sunspots  to  rotate  in  the 
same  direction  as  compared  to  opposite  directions  (Yan  et  al. 
2008). 

Multiple  causes  for  these  rotations  have  been  proposed,  but 
a  community  consensus  on  the  most  dominant,  or  sole,  mecha¬ 
nism  is  still  lacking.  Brown  et  al.  (2003)  suggest  photospheric 
flows  and  magnetic  flux  tube  emergence  as  potential  mecha¬ 
nisms,  while  Yan  &  Qu  (2007)  suggest  a  combination  of  the 
Coriolis  force  and  differential  rotation  as  alternative  mecha¬ 
nisms.  Su  et  al.  (2008)  investigated  the  interactions  of  subsur¬ 
face  vortical  flows  and  magnetic  flux  tubes  as  possible  drivers 
of  sunspot  rotation. 

Investigating  the  chromospheric  surface  above  sunspots, 
Balasubramaniam  et  al.  (2004)  have  shown  that  the  sunspot 
magnetic  structure  is  preserved  above  the  photospheric  layer, 
hence,  rotations  should  be  well  correlated  between  photosphere 
and  chromosphere. 

Brown  et  al.  (2003)  suggest  that  measured  sunspot  rotation 
rates  are  orders  of  magnitude  greater  than  what  would  be 
expected  from  rotations  caused  by  differential  rotation  or 
the  Coriolis  force.  Table  1  shows  a  sampling  of  the  average 
or  maximum  intrinsic  sunspot  rotation  rates  reported  in  the 
literature.  Note  that  most  reports  of  intrinsic  sunspot  rotations 
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Table  1 

Representative  Intrinsic  Sunspot  Rotation  Rates  Reported  in  the  Literature 


NOAA  Active  Region 

Peak  or  Average  Rotation  Rate 
(“hr-1) 

Peak  or  Average  Rotational  Velocity 
(km  s-1) 

Rotational  Direction 

Reference 

10759 

0.85  (avg.) 

ccw 

Kazachenko  et  al.  (2009) 

10930 

8  (peak) 

ccw 

Min  &  Chae  (2009) 

10930 

20  (peak) 

ccw 

Zhang  et  al.  (2007) 

8668 

1.3  (peak) 

0.055 

ccw 

Brown  et  al.  (2003) 

9004 

3  (peak) 

0.095 

cw 

Brown  et  al.  (2003) 

9077 

1.2  (peak) 

0.030 

ccw 

Brown  et  al.  (2003) 

9114 

2.2  (peak) 

0.077 

ccw 

Brown  et  al.  (2003) 

9280 

0.8  (peak) 

0.031 

ccw 

Brown  et  al.  (2003) 

9354 

1 .4  (peak) 

0.030 

ccw 

Brown  et  al.  (2003) 

10030 

2  (peak) 

0.077 

ccw 

Brown  et  al.  (2003) 

10424 

2.7  (avg.) 

ccw 

Yan  &  Qu  (2007) 

8210 

Few  degrees  per  hour 

cw 

Regnier  &  Canfield  (2006) 

10486 

1.2  (leading  spot  avg.) 

ccw 

Zhang  et  al.  (2008) 

10486 

1.5  (following  spot  avg.) 

ccw 

Zhang  et  al.  (2008) 

10486 

1.5  (leading  spot  avg.) 

0.200  (+0.023/— 0.028) 

cw 

This  work 

10486 

0.76  (following  spot  avg.) 

0.131  (+0.011/— 0.015) 

ccw 

This  work 

Notes.  Average  or  maximum  rotation  rates  are  reported.  Note  that  rotation  rates  are  usually  variable  through  time  when  rotation  occurs.  Rotational  velocities  are 
calculated,  when  possible,  using  Equation  (1)  from  Brown  et  al.  (2003).  For  the  leading  sunspot  in  NOAA  AR  10486,  Zhang  et  al.  (2008)  report  counter-clockwise 
rotation  on  five  of  the  six  days  the  sunspot  was  measured.  Small  clockwise  rotation  was  measured  on  2003  October  27.  On  2003  October  29,  Zhang  et  al.  (2008)  report 
that  the  leading  sunspot  rotated  at  an  average  angular  speed  of  4?  3  8  hr-1 . 


display  variable  rotation  rates  through  time  when  a  sunspot  is 
rotating  (Brown  et  al.  2003;  Min  &  Chae  2009). 

Another  question  involves  the  relationship  between  intrinsic 
sunspot  rotations  and  the  occurrence  of  flares.  If  a  causative 
connection  exists,  then  it  may  be  possible  to  predict  flaring, 
which  would  have  implications  for  space  weather  prediction 
efforts.  Stenflo  (1969)  and  Barnes  &  Sturrock  (1972)  suggest 
that  intrinsic  sunspot  rotation  can  lead  to  magnetic  energy 
build-up  and  release,  directly  followed  by  a  solar  flare.  Brown 
et  al.  (2003)  studied  seven  active  regions  and  compared  time- 
coincident  occurrences  of  sunspot  rotations  and  flares.  While 
flares  occurred  in  six  of  the  seven  active  regions  during  sunspot 
rotations,  Brown  et  al.  (2003)  suggest  that  only  three  of  the  flares 
are  likely  to  be  associated  with  sunspot  rotations.  Min  &  Chae 
(2009)  report  significant  rotation  of  the  small  sunspot  in  NOAA 
AR  10930,  but  no  flaring  has  been  causatively  linked  with  this 
rotation. 

Yan  &  Qu  (2007),  Zhang  et  al.  (2007),  and  Yan  et  al.  (2009) 
suggest  direct  relationships  between  rotating  sunspots  and  flares, 
while  Kazachenko  et  al.  (2009)  suggest  that  sunspot  rotation 
injected  significant  helicity  into  the  corona  of  NOAA  AR  10759, 
that  led  to  flaring.  Kazachenko  et  al.  (2010),  however,  report 
that  sunspot  rotation  in  NOAA  AR  10486  did  not  contribute 
significantly  to  the  helicity  and  energy  budgets  that  led  to 
some  of  the  major  solar  flares  during  the  2003  “Halloween 
storms.” 

The  summation  of  past  research  on  connections  between  in¬ 
trinsic  sunspot  rotation  and  flaring  requires  further  demonstra¬ 
tion  and  a  larger,  detailed  effort  to  search  for  systematics  in 
sunspot  rotational  behavior. 

Chromospheric  filaments  are  another  atmospheric  phenom¬ 
ena  that  have  associations,  and  possible  implications,  to  so¬ 
lar  flares.  Filament  mass  motion  and  flare  initiation  have  been 
studied  by  many  workers  (Jing  et  al.  2003,  2004;  Sterling  & 
Moore  2005;  Rubio  da  Costa  et  al.  2012;  Vemareddy  et  al.  2012; 
Zuccarello  et  al.  2009).  This  research  will  add  to  the  growing 
body  of  knowledge  concerning  both  intrinsic  sunspot  rotation 


and  filament  mass  motion,  which  often  contribute  to,  or  are 
associated  with,  solar  flares. 

2.  OBSERVATIONS  AND  METHODOLOGY 

Continuum  (630.32  nm)  and  Ha  (656.28  nm)  imagery  for 
this  research  were  obtained  from  the  USAF  ISOON  on  2003 
October  27,  28,  and  29.  At  the  time  of  the  data  set,  ISOON  was 
located  at  the  National  Solar  Observatory,  Sacramento  Peak, 
New  Mexico,  at  an  altitude  of  2939  m.  Primary  ISOON  op¬ 
tical  and  mechanical  components  include  a  250  mm  diameter 
objective,  independent  right  ascension  and  declination  motor 
drives,  tunable  0.1  A  Fabry-Perot  etalons,  and  Cooke  triplet 
lenses.  CCD  images  were  obtained  by  a  cooled  (  —  10°C)  2048  x 
2048  pixel  array  with  a  plate  scale  of  1.08437  arcsec  pixel-1. 
Original  ISOON  solar  images  have  a  plate  scale  of 
781  km  pixel-1  or  720  km  arcsec-1 .  The  plate  scale  for  flattened 
images  is  498  km  pixel-1  or  459  km  arcsec-1 .  Dark  current  and 
flat-field  calibration  frames  were  automatically  applied  to  all  im¬ 
ages  immediately  following  image  acquisition.  Figure  1  shows 
an  example  of  original  and  flattened/limb-darkening-removed 
ISOON  continuum  and  Ha  images. 

Continuum  images  on  October  27  and  29  had  a  cadence  of 
30  minutes  or  1  hr  ranging  from  ~  14:00  to  23:30  UT  or  within  a 
subset  of  that  timeframe.  Sets  of  three  continuum  images  were 
taken  during  each  active  observation  with  the  highest-resolution 
image  used  for  analysis.  Only  five  sets  of  continuum  images  on 
October  28  were  available,  which  did  not  allow  analysis  of 
photospheric  motions  for  that  day.  Ha  images  with  a  1  minute 
cadence  were  available  from  ~  14:00  to  23:00  UT  on  all  three 
days.  Data  gaps  on  October  29  from  17:22  to  18:01  UT,  19:46 
to  20:03  UT,  20:13  to  21:01  UT,  and  23:30  to  23:46  UT  exist  in 
the  Ha  data  set  due  to  clouds.  Ha  image  analysis  used  images 
with  a  10  minute  cadence. 

NOAA  Active  Regions  (AR)  10484,  10486,  and  10488  were 
the  most  prominent  active  regions  visible  in  the  data  set  and  are 
the  focus  of  this  paper.  Data  analysis  techniques  applied  to  these 
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Figure  1.  Representative  ISOON  continuum  images  (top  row)  and  Ha  images  (bottom  row)  from  2003  October  27.  The  images  on  the  left  are  the  original  ISOON 
images.  The  images  to  the  right  represent  the  images  on  the  left  after  being  flattened  with  solar  limb  darkening  removed. 


active  regions  include  sunspot  uncurling  to  measure  rotational 
motions  in  the  photosphere  and  chromosphere,  and  detection  of 
filament  motion  during  active  region  evolution  via  IDL  movie 
routines  and  SAOImage  DS9  (Joye  &  Mandel  2003;  Joye  2006). 

Uncurling  sunspots  to  measure  rotational  motion  follows  the 
techniques  described  in  Brown  et  al.  (2003).  Using  IDL  routines, 
the  ISOON  images  chosen  for  analysis  were  first  flattened  with 
limb  darkening  removed.  In  each  set  of  continuum  and  H a 
images,  each  active  region  was  extracted  by  defining  a  large 
rectangular  area  around  the  active  region.  The  cropped  set  of 
images  for  each  active  region  were  then  aligned  relative  to 
a  reference  image,  normalized,  and  trimmed  to  remove  edge 
artifacts. 

Both  large  sunspots  in  NOAA  AR  10486  and  10488  were 
uncurled,  while  the  leading  and  trailing  sunspots  of  the  four- 
sunspot  complex  in  NOAA  AR  10484  were  uncurled.  For  each 
sunspot,  a  circular  aperture  was  centered  on  each  sunspot  of 
a  chosen  radius  that  extends  well  beyond  the  penumbra.  Each 
sunspot  was  then  uncurled — transformed  from  an  (x,  y)  to  (r,  0) 
coordinate  system — in  1°  increments  beginning  with  the  eastern 
chord  (0°)  and  moving  sequentially  CCW.  Figure  2  shows  an 
uncurled  image  for  the  following  sunspot  of  NOAA  AR  10486. 

Measuring  the  rotational  and  dynamical  motion  of  sunspots 
and  filaments  requires  extracting  radial  time  slices  from  un¬ 
curled  images  and  measuring  non-horizontal  motion  in  ( 0 ,  t) 


Figure  2.  A  representative  uncurled  continuum  image  ( 6 ,  r)  of  the  following 
large  sunspot  in  NOAA  AR  10486  on  October  29.  The  image  shows  the  dark 
umbra  stretched  across  the  first  ~15  radial  pixels  and  the  asymmetrical  leading 
sunspot  at  ~  1 80°  along  with  its  surrounding  penumbra  and  pores.  Some  artifacts 
from  the  reduction  process  are  present  along  the  y-axis.  See  the  text  for  details 
of  the  image  uncurling  process. 
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Figure  3.  (a)  An  angle  vs.  time  plot  for  the  following  sunspot  in  NOAA  AR  10486  on  October  29.  Time  slices  10  pixels  wide  were  extracted  from  continuum  images 
and  aligned  chronologically.  Note  the  wide  area  of  rotational  (i.e.,  non-horizontal)  motion  from  ~230°  to  300°.  The  time  range  for  the  plot  is  from  ~14:15UT 
to  ~22:30UT.  (b)  Plot  showing  the  measurement  of  two  of  the  features  within  the  area  of  rotation  in  panel  (a).  Error  bars  represent  standard  deviation  from 
10  measurements  of  the  center  of  each  feature. 


space  (see  Figure  7;  Brown  et  al.  2003).  Features  that  exhibit 
a  change  in  positional  angle  through  time  (A0 /At)  exhibit  rota¬ 
tional  motion,  while  those  not  showing  such  angular  movement 
do  not.  Four  radial  slices  with  10  pixel  widths  (10/84),  rang¬ 
ing  from  the  inner  penumbra  to  the  outer  penumbra/quiet  pho¬ 
tosphere,  were  extracted  from  each  uncurled  sunspot,  aligned 
chronologically  on  each  day,  and  examined  for  rotational  mo¬ 
tion  (in  ( 0 ,  t)  space).  The  pixel  coordinates  for  each  feature 
(continuum  and  Ha)  were  manually  measured  across  the  times- 
pan  of  the  data  to  determine  its  position,  0 ,  which  allowed  the 
calculation  of  the  feature  rotation  rate  in  0  hr-1  and  km  s-1  via 
Equation  (1)  in  Brown  et  al.  (2003).  Ten  measurements  were 
made  for  each  time  slice  with  the  mean  and  standard  deviation 
calculated  to  determine  feature  position  through  time.  Figure  3 
shows  both  a  time- slice  plot  of  the  following  sunspot  rotation 
on  October  29  for  NOAA  AR  10486  and  the  two  features  in  that 
region  that  were  measured. 

IDL  movie  routines  and  SAOImage  DS9  were  also  used  to 
study  time-lapse  movies  of  the  three  active  regions.  The  primary 
benefit  of  this  effort  was  to  correlate  chromospheric  filament 
motions  in  each  active  region  with  the  rotational  motions  seen 
in  the  Ha  time- slice  plots  of  the  uncurled  sunspots.  Movies  were 
also  beneficial  in  detecting  very  fast  and/or  complex  filament 
motion  that  is  not  detectable  in  the  time-slice  plots. 

Positive  and  negative  magnetic  field  areas  in  SOHO  magne¬ 
togram  images  were  measured  through  time  using  the  IDL  pro¬ 
gram,  Helioflux  (Panasenco  et  al.  2011).  This  was  conducted 
to  understand  the  broader  magnetic  field  context  in  each  active 
region  through  the  time  of  the  data  set. 


3.  RESULTS 
3.1.  NOAA  AR  10484 

NOAA  AR  10484  became  visible  on  the  solar  surface  on 
October  18  and  remained  visible  until  rotating  out  of  view 
on  October  31.  At  its  peak,  this  AR  had  a  magnetic 
classification  (http://solar.ifa.hawaii.edu/ARMaps).  Continuum 
ISOON  images  from  October  27  reveal  an  active  region  with  five 
or  six  discernible  umbrae  of  various  sizes  with  penumbra  linking 
the  umbrae  together.  Ha  images  on  October  27  clearly  show 
the  sunspots  surrounded  by  diffuse  chromospheric  brightening 
around  the  active  region,  dominantly  to  the  southwest  of 
the  sunspots.  Visible  filaments  are  present  to  the  northeast, 
southeast,  south,  and  northwest  of  the  sunspots. 

Two  C-class  solar  flares  occurred  in  our  data  set  near 
AR  10484  on  October  27,  including  a  C7.5  flare  that  be¬ 
gan  at  18:34  UT  and  a  C9.0  flare  that  began  at  19:48  UT 
(http://www.lmsal.com/solarsoft/last_events_20031028_0955/ 
index.html).  The  C7.5  flare  shows  Ha  brightening  to  the  south 
of  the  leading  sunspot  and  to  the  northeast  of  the  sunspots, 
while  Ha  brightening  for  the  C9.0  flare  occurred  solely  to  the 
northeast  of  the  active  region. 

There  was  no  rotational  motion  apparent  in  the  continuum 
images  on  October  27  for  either  the  leading  or  following  sunspot. 
Two  dark  chromospheric  features  in  time-slice  sequences  at 
average  distances  of  45  pixels  (or  48/8)  from  the  umbral  center 
were  seen  in  uncurled  images  of  the  following  sunspot  at  ~  145° 
and  ~290°,  respectively.  The  feature  at  ~  145°  displayed  motion 
from  ~16:33UT  to  ~21:20UT,  while  the  feature  at  ~290° 
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17:43:04  UT 


Figure  4.  Time  sequence  of  four  Ha  images  of  NOAA  AR  10484  on  2003 
October  27,  from  16:13  UT  to  17:43  UT.  Discrete  filament  mass  can  be  seen 
migrating  north  in  a  reverse  sigmoid  fashion,  beginning  with  the  top  image  and 
moving  downward.  Additional  mass  movement  at  ~  19:30  UT  directly  precedes 
a  C9.0  class  flare  that  begins  at  ~19:44UT. 

(An  animation  of  this  figure  is  available  in  the  online  journal.) 


showed  rotational  motion  from  ~16:33  UT  to  ~19:45  UT.  Both 
features,  which  are  likely  to  be  associated  with  filaments,  moved 
CCW  at  average  rates  of  2?28hr-1  (0.39  km  s-1)  and  UlOhr-1 
(0.19kms_1). 

Time-series  analysis  of  AR  10484  Ha  images  in  SAOImage 
DS9  (Joye  &  Mandel  2003;  Joye  2006)  shows  a  series  of  discrete 
filament  motions  north/northeast  of  the  following  sunspot  prior 
to  both  the  C7.5  flare  and  the  C9.0  flare.  Filament  material 
followed  a  reverse- S  path  that  begins  just  north  of  the  following 
sunspot  at  —  16: 13  UT  and  moves  in  a  distinct  reverse-S  pattern 
north  of  the  same  sunspot  by  —  17:43  UT.  Figure  4  displays  a 
sequence  of  images  tracing  the  path  of  the  filamentary  mass, 
which  can  be  seen  traveling  up  to  —95"  (—44,000  km)  from 
the  following  sunspot  (a  movie  of  this  motion  is  provided  with 
the  online  version  of  this  paper).  This  movement  precedes  the 


occurrence  of  the  C7.5  flare  by  — 2  hr  and  is  spatially  coincident 
with  brightening  resulting  from  the  flare. 

A  second  burst  of  filament-like  mass  motion  occurred  from 
— 19:14  UT  to  — 19:24  UT  in  a  direction  that  is  halfway  along 
the  reverse-S  path  moving  to  the  north  (about  the  same  position 
as  the  17:13  UT  image  in  Figure  4).  The  path  of  both  sets  of 
filamentary  motion  is  spatially  coincident  and  the  second  set 
of  mass  motion  just  precedes  the  C9.0  flare,  where  brightening 
occurs  at  ~  19:48  UT.  Filamentary  mass  motion  occurred  along 
the  northwest  boundary  of  the  Ha  brightening  from  the  C9.0 
flare  prior  to  the  occurrence  of  the  flare. 

On  October  28,  there  was  no  detectable  photospheric 
mass  motion  in  continuum  images.  Ha  chromospheric  images 
showed  filamentary  material  moving  CCW  at  an  average  ro¬ 
tational  speed  of  12?8hr_1  (1.17kms_1)  from  — 14:30  UT  to 
— 22:30  UT.  The  rotational  rates  were  variable  with  somewhat 
faster  rotation  occurring  from  — 14:30  UT  to  — 18:40  UT  and  a 
slower  rotational  rate  from  — 18:40  UT  to  — 22:30  UT.  Filament 
motion  on  October  28  must  have  occurred  by  a  mechanism  other 
than  sunspot  rotation  due  to  the  lack  of  any  detectable  rotation. 
Figure  5  displays  the  rotational  characteristics  of  the  filamentary 
motion  on  October  27  and  28. 

3.2.  NOAA  AR  10486 

NOAA  AR  10486  was  large  and  active  during  the  time 
range  of  the  data  set  and  one  of  the  principal  contributors  to 
the  “Halloween  solar  storms”  in  late  2003  October  and  early 
November  (Hady  2009).  AR  10486  emerged  on  the  solar  disk 
late  on  October  23  as  a  large,  mature  active  region.  A  large 
X-class  flare  occurred  during  the  three  days  of  our  data — an 
X10  flare  on  October  29  that  began  at  — 20:37  UT.  This  flare 
occurred  the  day  after  an  — X17  flare  that  began  at  ~9:5 1  UT 
on  October  28  that  has  been  extensively  studied  (Ambastha 
2007;  Zuccarello  et  al.  2009;  Kazachenko  et  al.  2010),  but  is 
not  available  in  our  data.  SOHO  magnetogram  measurements 
of  AR  10486  for  the  three  days  show  increasing  and  decreasing 
negative  and  positive  magnetic  field  areas,  respectively,  prior 
to  October  27,  and  the  leveling  off  of  the  magnetic  field  areas 
early  on  October  28.  Magnetic  field  areas  of  both  polarities 
remain  constant  through  the  rest  of  the  data  set  with  the  negative 
magnetic  field  area  —  1 .3  times  larger  than  the  positive  magnetic 
field  area. 

On  October  29,  significant  photospheric  and  chromospheric 
mass  motion  occurs  before,  during,  and  after  the  XI 0  flare 
that  is  prominent  in  the  Ha  images.  This  motion  is  mostly 
consistent  with,  and  a  subset  of,  the  five  days  of  rotational 
motion  seen  in  NOAA  AR  10486,  as  reported  by  Zhang  et  al. 
(2008).  Continuum  images  of  the  following  sunspot  show  a  wide 
swath  of  CCW  photospheric  mass  motion  from  -220°  to  -300° 
(Figure  3(a)).  Three  distinct  features  within  the  overall  area  of 
motion  were  measured  at  an  average  distance  of  48"8  (45  pixels) 
from  the  umbral  center  and  show  mass  motion  with  rotational 
velocities  of  0?70hr-1  (0.12kms_1),  0?71hr_1  (0.12kms_1), 
and  0?86hr-1  (0.15kms_1)  (Figure  3(b)).  Using  an  average 
photospheric  density  of  10-9  g  cm-3,  a  vertical  photosphere 
thickness  of  500  km  (Eddy  1979),  and  gas  moving  across  -70°, 
the  approximate  mass  of  photospheric  gas  entrained  in  rotation 
from  -14:20  UT  to  -22:30  UT  was  -7  x  1017  kg. 

In  contrast,  continuum  images  of  the  leading  sunspot  show 
a  much  more  restricted  region  of  CW  mass  motion  occurring 
from  — 14:00  UT  to  —  21:00  UT  with  an  average  velocity  of 
1?50  hr-1  (0.20  km  s-1).  The  feature  was  measured  at  an  average 
distance  of  38"  (35  pixels)  from  the  umbral  center.  This  feature 
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Figure  5.  Rotational  motion  of  filaments  as  measured  from  time  slices  of  the  following  sunspot  in  NOAA  AR  10484  on  2003  October  27,  top,  and  on  October  28, 
bottom.  The  motion  on  2003  October  27  corresponds  with  filament  mass  motion  moving  to  the  north/ northeast  of  the  following  sunspot  that  is  time  coincident  with  a 
C9. 0-class  solar  flare.  Vertical  bars  in  the  top  figure  represent  the  approximate  times  of  the  C7.5  and  C9.0  flare,  respectively. 
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has  been  identified  as  a  distended  segment  of  the  leading  large 
sunspot’s  umbra,  which  extends  to  the  west/ southwest  of  the 
primary  umbra.  This  motion  is  in  the  opposite  direction  to  that 
reported  by  Zhang  et  al.  (2008),  but  involves  different  areas  of 
the  sunspot. 

Ha  images  of  the  following  sunspot  show  systematic 
CCW  motions  (through  — 3°-4°)  of  bright  regions  in  the 
—14:20-17:00 UT  time  frame  at  average  rates  of  — 0?67  hr-1 
(0.14kms-1)  and  —  l?55hr-1  (0.33 kms-1).  Rotation  of  both 
features  ceases  after  — 18:00  UT  until  the  onset  of  the  X 10  flare 
at  — 20:37  UT.  Two  features  in  the  region  of  flare  brightening 
of  the  following  sunspot  appear  to  shift  abruptly  CCW  by  -8° 
and  -14°  before  returning  to  relatively  constant,  but  displaced, 
positions.  Ha  images  of  the  leading  sunspot  also  show  a  CCW 
shift  of  material  by  -10°  before  moving  CW  after  the  flare  to 
return  to  within  -2°  of  the  pre-flare  position. 

3.3.  NOAAAR  10488 

NOAA  AR  10488  appears  in  SOHO  magnetograms  on 
October  26  and  grew  rapidly  in  the  time  period  of  the  data  set. 
By  October  27,  continuum  images  revealed  an  elliptical  region 
outlining  development  of  the  leading  and  trailing  sunspots,  with 
the  leading  sunspot  incompletely  formed  with  only  a  partial 
umbra  and  penumbra.  Ha  images  showed  the  active  region 
enveloped  in  arc  filaments  stretching  from  the  leading  sunspot 
to  the  following  sunspot. 

The  large,  leading  sunspot  was  fully  formed  by  October  28  in 
continuum  images,  with  the  following  spot  continuing  to  grow 
via  the  accumulation  of  pores  throughout  October  28  and  29. 
The  leading  sunspot  also  became  visible  in  Ha  images  on 
October  28,  while  the  following  sunspot  was  only  partially 
visible  on  October  28  and  29.  Positive  and  negative  magnetic 
field  areas  grew  at  similar  rates  during  the  three  days,  with  the 
magnetic  field  areas  increasing  in  size  by  an  order  of  magnitude 
(— 109  to  — 1010  km2)  between  October  26  and  31. 

Arc  filaments  for  both  sunspots  on  October  27  show  both 
CW  and  CCW  motion  at  varying  rates  and  locations  around 
the  sunspots,  and  over  different  time  intervals  in  Ha  images. 
Two  filamentary  features  around  the  leading  sunspot  show 
fairly  regular  CCW  motion  from  — 14:15  UT  to  —  16:50  UT  at 
rates  of  5?23hr-1  (0.90 kms-1)  and  l?46hr-1  (0.31  kms-1). 
Another  filamentary  feature,  at  -315°,  moved  CW  at  1?43  hr-1 
(0.30  kms-1). 

Filamentary  mass  motion  around  the  following  sunspot  on 
October  27  also  showed  CW  and  CCW  motions.  Two  fila¬ 
mentary  fragments,  measured  at  different  radial  distances  from 
the  umbral  center,  likely  represent  the  same  filament.  One  seg¬ 
ment,  48.8”  from  the  umbral  center,  moved  abruptly  CW  from 
-20:45 UT  to  -22:00 UT  at  3?43hr-1  (0.59 kms-1),  while 
the  segment  measured  at  59.6”  showed  irregular  CW  motion 
from  — 18:15  UT  to  — 22:00 UT  at  average  rates  of  l?63hr-1 
(0.34  km  s-1).  A  filament  at  a  different  location  around  the  fol¬ 
lowing  sunspot  moved  consistently  CCW  from  — 15:30  UT  to 
— 23:00  UT  at  2?  22  hr-1. 

Filaments  surrounding  the  leading  sunspot  dominate  the  mo¬ 
tion  seen  on  October  28.  One  filament,  moving  at  a  CCW  rate  of 
4?35  hr-1  (0.91  km  s-1),  only  shows  motion  from  —14:30  UT  to 
— 16:30  UT.  Three  different  filaments,  however,  show  move¬ 
ment  from  — 14:30  UT  to  — 22:30 UT  at  average  rates  of 
0?66 hr-1  (0.11  kms-1:  CCW),  0?97hr-1  (0.17 kms-1:  CW), 
and  1?23 hr-1  (0.21  kms-1:  CW). 

Two  individual  pores  trace  CW  and  CCW  rotational  motion 
in  continuum  images  on  October  29  as  they  move  away  from 
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Figure  6.  Time-slice  plot  of  the  open  photosphere  to  the  west  of  the  leading 
sunspot  in  NOAA  AR  10488  on  October  29,  2003.  The  time  span  of  the  data 
ranges  from  ~14: 14  UT  to  ~22:30  UT.  Two  pores  are  approaching  each  other  at 
~180°  and  appear  to  be  in  the  process  of  merging.  The  black  horizontal  dashed 
lines  are  artifacts  of  the  data  reduction  process. 


the  leading  sunspot.  Figure  6  shows  a  time- slice  plot  of  the  two 
pores  at  -170°  where  the  top  pore  is  moving  CW  at  an  angular 
rate  of  0?57hr-1  (0.12  km  s-1)  while  the  bottom  pore  moves 
CCW  at  a  rate  of  0?38  hr-1  (0.07  km  s-1).  Two  pores  around  the 
following  sunspot  were  moving  CW  at  rates  of  2?  15-2?  31  hr-1 
(0.3 1-0.37 kms-1)  while  a  third  pore  was  moving  CCW  at 
0?57hr-1  (0.12  kms-1). 

Measurable  filamentary  motion  on  October  29  was  restricted 
to  the  following  sunspot,  although  some  short-duration  motion 
to  the  southwest  of  the  leading  sunspot  was  seen.  Filamentary 
motion  occurred  primarily  to  the  southeast  of  the  following 
sunspot  and  was  visible  throughout  the  day’s  images.  One  fil¬ 
ament,  measured  at  average  distances  of  38"8  and  48"8  from 
the  umbral  center,  displays  motion  at  rates  of  1?  16-2?  16  hr-1 
(0.16-0.37 kms-1:  CCW).  Another  filament,  -60°  away, 
rotated  CW  at  l?06hr-1  (0.22 kms-1). 

The  photospheric  pore  motions  and  the  filament  motions 
measured  in  this  young  active  region  are  not  related,  but  are 
typical  of  both  photospheric  and  chromospheric  activity  in  a 
growing  bi-polar  sunspot  pair  and  active  region. 

4.  ANALYSIS  AND  CONCLUSIONS 

From  the  work  of  Zhang  et  al.  (2008)  and  this  work,  it  appears 
that  countervailing  rotational  motions  are  possible  around  the 
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Figure  7.  Episodic  CCW  motion  of  chromospheric  mass  of  two  features  in  the  following  sunspot  of  NOAA  AR  10486  during  and  immediately  following  an  X10  flare, 
which  is  indicated  by  the  offset  of  two  chromospheric  features  on  2003  October  29  at  ~20:45  UT.  The  X10  flare  began  at  20:37  UT  and  ended  at  21:01  UT,  which  is 
indicated  by  the  vertical  lines  in  both  figures.  The  location  of  this  mass  movement  is  near,  or  adjacent  to,  the  mass  motion  in  the  photosphere. 


same  sunspot  at  the  same  time.  Zhang  et  al.  (2008)  report  CCW 
rotation  on  October  29  for  both  sunspots  in  NOAA  AR  10486, 
while  this  work  identified  an  extended  umbral  segment  to  the 
west/ southwest  of  the  leading  sunspot,  embedded  in  penumbra, 


that  exhibits  CW  rotation.  The  small  region  of  umbral  CW 
motion  seen  in  this  work  has  a  negative  magnetic  polarity  that 
directly  abuts  the  region  of  large  positive  polarity  for  the  leading 
sunspot,  as  seen  in  SOHO  magnetograms. 
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Another  interesting  issue  is  the  lack  of  large-scale,  complete 
rotational  penumbral  motion  around  a  sunspot  umbra.  Previous 
works  have  identified  various  small  features  around  a  sunspot 
that  exhibit  rotation  (Brown  et  al.  2003;  Yan  &  Qu  2007;  Min 
&  Chae  2009;  Kazachenko  et  al.  2009,  2010;  Yan  et  al.  2009; 
Zhang  et  al.  2007).  However,  there  has  yet  to  be  a  report  of  full 
penumbral  rotation  around  a  sunspot  center.  The  extent  of  mass 
entrainment  seen  in  the  following  sunspot  in  NOAA  AR  10486 
on  October  29  is  significant  and  hints  at  a  region  with  magnetic 
fields  of  sufficient  strength  to  entrain  the  surrounding  plasma 
in  the  photosphere.  The  lack  of  homogeneous  360°  rotation 
around  a  sunspot  suggests  variable  magnetic  field  strengths, 
which  may  be  detectable  and  measurable  in  extreme  ultraviolet/ 
X-ray  coronal  images. 

The  region  of  mass  motion  in  the  following  sunspot  in 
NOAA  AR  10486  on  October  29  displays  linear  to  semi-linear 
motion  during,  and  after,  the  X10  flare.  No  major  rotational 
acceleration  or  deceleration  resulted  from  the  occurrence  of 
the  flare.  Kazachenko  et  al.  (2010)  suggests  that  rotation  in 
AR  10486  on  October  28  provided  a  negligible  component  to  the 
helicity  and  energy  budget  in  the  region  where  flaring  occurred. 
If  this  condition  persisted  into  October  29,  then  that  may  be  why 
there  was  little  or  no  rotational  response  during  or  after  the  X10 
flare. 

Lastly,  NOAA  AR  10486  displays  some  indications  of 
the  mass  coupling  of  the  photosphere  and  chromosphere  on 
October  29,  primarily  before  the  X10  flare.  While  the  rotation 
of  regions  of  the  photosphere  is  prominent  (Figure  3)  in  the  fol¬ 
lowing  sunspot,  CCW  chromospheric  motion  is  much  smaller 
in  magnitude  and  more  episodic.  Some  chromospheric  mass 
also  appears  to  have  been  displaced  as  a  result  of  the  X10  flare 
(Figure  7).  This  chromospheric  motion  is  time  coincident  with 
the  resulting  proton  storm  as  seen  in  SOHO  EUV  Imaging  Tele¬ 
scope  images  of  the  same  time  period.  The  abrupt  filament 
movements  in  close  association  with  the  XI 0  flare  could  have 
been  caused  by  a  Ha  Moreton  wave  or  a  magnetosonic  shock 
wave  propagating  through  the  lower  regions  of  the  solar  atmo¬ 
sphere  (Shen  &  Liu  2012;  Asai  et  al.  2012). 

One  suggestion  for  the  observed  photospheric- 
chromospheric  coupling  is  the  presence  of  magnetic  field  lines 
of  varying  inclinations  from  the  vertical  that  are  preferentially 
stronger  in  the  penumbral  regions  of  mass  motion  compared  to 
the  field  lines  originating  from  other  locations  in  the  penum¬ 
bra.  While  it  is  well  known  that  sunspot  magnetic  field  strength 
decreases  with  increasing  radial  distance  from  the  sunspot  cen¬ 
ter  (Bhatnagar  &  Livingston  2005),  this  result  suggests  greater 
variations  in  magnetic  field  strength  as  a  function  of  position 
around  the  entire  penumbra. 

Ha  images  in  this  work  sample  the  core  of  the  feature  at 
656.28  nm  (with  a  bandwidth  of  —0.080  A).  Radiative  transfer 
and  model  atmosphere  calculations  for  the  disk  chromosphere 
show  that  the  core  of  the  Ha  line  samples  the  lower  to 
middle  portions  of  the  average  chromosphere  (1600-3000  km: 
White  &  Wilson  1966;  1200-1700  km:  Vernazza  et  al.  1981; 
1107-2100  km:  Qu  &  Xu  2002;  1900  km:  Balasubramaniam 
et  al.  2007;  500-1500  km:  Leenaarts  et  al.  2012).  Hence, 
we  estimate  that  the  measurements  here  represent  the  average 
photosphere  and  chromosphere  and  are,  at  best,  separated  by  a 
few  thousand  kilometers. 

We  find  it  important  to  address  the  possible  application  of  non¬ 
linear  force-free  (NLFF)  modeling  methods  to  understand  the 
rotational  coupling  in  the  photosphere  and  chromosphere.  NLFF 
models  (see,  e.g.,  Wiegelmann  &  Sakurai  2012;  Wiegelmann 
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et  al.  2012)  are  still  in  their  infancy.  The  vertical  resolution/ 
granularity  of  the  extrapolation  is  not  currently  estimated 
(perhaps  a  few  thousand  kilometers)  and  never  verified.  It  is 
not  a  self-consistent  problem  where  one  can  fit  the  observed 
magnetic  fields  to  vertical  structures. 

Modeling  of  the  magnetic  field  requires  vector  magnetic 
measures  that  are  limited  to  one  layer,  namely  the  photosphere, 
which  is  also  its  lower  bound.  We  are  still  unsure  how  surface 
errors  of  the  measured  longitudinal  field  (up  to  20%  for  weak 
fields,  e.g.,  a  5  G  error  in  a  20  G  measurement)  and  transverse 
fields  (—50%,  e.g.,  a  100  G  error  in  a  200  G  measurement) 
translate  to  an  error  in  the  vertical  direction  as  vertical  fields  are 
assumed  to  have  an  exponential  drop  in  the  modeling  processes. 
Measuring  the  change  in  the  magnetic  field  with  height  in  the 
solar  atmosphere  (i.e.,  dB/dz )  also  requires  accurate  measures 
of  gas  pressure  and  magnetic  pressure  variations  with  altitude. 
Consistent  measurements  for  photospheric  magnetic  fields  are 
now  common  (Petrie  2012),  but  chromospheric  magnetic  field 
measurements  are  still  lacking.  Hence,  using  extrapolation 
methods  is  unproven  to  explain  the  twisting  of  sunspots  in  two 
layers. 

We  plan  to  further  analyze  similar  data  sets  to  understand  any 
coupling  relationships  between  the  photosphere  and  chromo¬ 
sphere,  using  both  ground-based  ISOON  and  Solar  Dynamics 
Observatory  data. 

We  are  grateful  for  the  use  of  the  ISOON  data.  P.S.H.  and  S.S. 
were  supported  by  UND  for  this  work.  K.S.B.  was  supported  by 
the  AFOSR  and  AFRL  in  this  work. 
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